Manipulating the fully spin-polarized edge currents in graphene ribbon 
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Electron fully spin-polarized edge states in graphene emerged at the interfaces of a nonuniform 
magnetic field are studied numerically in a tight-binding model, with both the orbital and Zeeman- 
splitting effects of magnetic field considered. We show that the fully spin-polarized currents can 
be manipulated by a gate voltage. In order to make use of the fully spin-polarized currents in the 
spin related transport, a three-terminal experiment is designed and expected to export the fully 
spin-polarized currents. This may have important applications in spin based nanodevices. 
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I. INTRODUCTION 



Recently, spin-dependent transport phenomena have 
received intensive studies because of their potential ap- 
plications in spintronics devices. Understanding how to 
generate and manipulate the spin-polarized currents is 
one of the key points for the development of the spin- 
based devices. Rashba spin-orbit interaction, which can 
be controlled by an (or effective) electric field, provides an 
efficient way to achieve a spin-polarized current, for ex- 
ample in graphene ribbonji doped semiconductors^ and 
2D electron systemi^ In most of these applications, the 
induced charge currents are only partially polarized,— but 
a fully spin-polarized (FSP) current is much more diffi- 
cult to achieve experimentally. Up to now, FSP cur- 
rent is found to be realized only in very limited materi- 
als, for example, in the half-metallic states of the ferro- 
magnetic metals^i^ and ferromagnetic semiconductors 
where the electronic valance band is partially filled 
for one spin component, whereas completely empty for 
the other. From material viewpoint, two-dimensional 
(2D) graphene, candidate material of new genera- 
tion electronics^rJ^ for its high mobility and low carrier 
density, is expected to have great application in spin- 
tronics devicesjii"— How to realize a FSP current in 
graphene is then quite meaningful for future graphene- 
based electronic devices. In graphene nanoribbon with 
zigzag edges, it has been predicted from a fist-principle 
calculation^^ and numerical simulationsi^iii that with 
the application of a transverse electric field, a half- 
metallic state can be achieved, which leads naturally to 
FSP edge currents propagating along sample boundaries. 

On the other hand, for a 2D electron system, a suffi- 
cient strong magnetic field will polarize the spins of all 
electrons in the same direction as the field, but only a 
much smaller magnetic field is needed to achieve a FSP 
current*^ The cost is that a nonuniform magnetic field 
(NMF) must be applied to the system. Actually, it has 
been predicted that a partially spin-polarized current 
will appear at the interfaces of inhomogeneous magnetic 
fields i^^'^° but much less attention are taken on the FSP 
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FIG. 1: (Color online) (a) Schematic illustration of the sys- 
tem: an array of ferromagnetic strips with different magneti- 
zations is placed on the top of a rectangular graphene layer. 
In the two-strip case, (b) the magnetic field (solid line) and its 
gauge potential (dotted line) as functions of x, and (c) the cor- 
responding energy spectrum for a zigzag-edge graphene rib- 
bon with L = 213 nm, / = 2.13 nm, Bi = 2 T and B2 = -3 T. 



states at the interfaces, especially, in graphene. 

In this paper, we study monolayer graphene in the 
presence of a stepwise NMF, which can be realized ex- 
perimentally in several ways?^"— to exhibit systemati- 
cally how to manipulate the FSP currents generated at 
the interfaces of the NMF. The situation is schematically 
shown in Fig. [ija). For the uniform magnetic field case, 
the FSP state is found to be absent, and the transport 
in graphene can be described by the so-called counter- 
propagating chiral edge states^^ when the chemical po- 
tential fi is between the spin gap due to Zeeman splitting. 
However, for the stepwise NMF case, several spin gaps 
are opened and the FSP edge states may emerge at the 
magnetic interfaces when ^ is between the gaps. In the 
following we will address the approach to manipulate this 
novel FSP edge currents. Furthermore, an experimental 
setup is designed to export the edge currents. 
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II. MODEL 



We begin our discussion by the tight-binding model on 
a honeycomb lattice in the presence of a perpendicular 
magnetic field B = (0, 0, Bz{x)), which is given by, 



<ij> 



(1) 

where the first term describes electrons hopping be- 
tween the nearest neighbors, suffering an additional 
phase caused by the orbital effect of magnetic field, with 
cj = {cl^jc]^) electron creation operators, whereas the 
second one the Zeeman splitting term. Here hz{x) = 
^g^sBzix) with g as the Lande factor. The magnetic 
field Bz {x) is assumed to change only in x direction, and 
the stepwise one which is focused on in this paper, has 
the following form. 



Bi < X < L, 

Bzix) = { Bi + ^2^{x - L) L<x<L + l, 

B2 L + l <x <2L + l. 

(2) 

For the gauge potential A, the Landau gauge A = 
{0,Ay{x), 0) is adopted, with Ay{x) = Bz{x')dx' given 
as: 



Bix 
Ay{x) = iBix 



< X < L, 
^(.t2-L2) L<x<L + l, 



B2X - ^^^^H L + l <x <2L + l. 

(3) 

For simplicity, in the following discussion, energy is mea- 
sured in unit of the hopping integral t. 
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FIG. 2: (Color online) Top panels: Electron energy spectrum 
of graphene ribbon under a stepwise magnetic field with one 
interface, with L = 213 nm (1000 zigzag chains) for (a) Case 
I with Bi = 2 T and B2 = 3 T, (b) Case II with Bi = 2 T 
and B2 = —3 T. Here the transition length for the varia- 
tion of magnetic field long x direction is I — 2.13 nm (10 
zigzag chains). The red and black lines denote bands of spin- 
up and spin-down states, respectively. The interface edge 
states are labelled with thick lines. The Zeeman energy is 
chosen to be one tenth of the largest LL spacing of the region 
with B — 2 T. The black dashed lines give the representative 
positions of the chemical potentials fi. Other panels: The 
corresponding four-terminal configurations in Hall measure- 
ment. The spatial distribution of the edge currents is shown, 
where the red and black lines with arrows represent the spin- 
up and spin-down electron edge currents, respectively. The 
thick lines with arrow at the interfaces in (b3) represent that 
each current is carried by two edge states. 



III. RESULTS AND DISCUSSIONS 

To analyze all possible edge states, we choose open 
(periodic) boundary condition in the x (y) direction, and 
numerically diagonalize the Hamiltonian Hq on a rectan- 
gular sample under a NMF to obtain the electron energy 
spectrum. In one interface case, the spectrum is shown in 
Fig-HJc), which is symmetric with respect to zero energy 
due to the particle-hole symmetry ( which, strictly speak- 
ing, should be combined with the time-reversal symme- 
try) preserved even in a NMF. The spectrum is composed 
of two scries of Landau levels (LLs) , including all the bulk 
states of the system, and branches of edge states, which 
play a dominant role in quantum transport in the pres- 
ence of magnetic field. The edge states can be classified 
into the normal ones located at the sample boundaries, 
and the interface ones located at the interface where the 
magnetic field changes abruptly. 



A. Generation and manipulation of the FSP edge 
currents 

Let's focus our discussion on the regime where the 
chemical potential fi is between the zcroth LLs of the 
two spin species. To make things more explicit, we show 
in Figs, mja) and [Hb) the energy spectrum around zero 
energy. The splitting of the zcroth LL and formation of 
two spin gaps Ai and A2 can be clearly seen. When fi 
lies within the interval |^| < Ai/2 or > A2/2, there 
exist only sample boundary edge states, and no interface 
edge states at all. However, when jj, lies between the two 
spin gaps Ai/2 < < A2/2, interface edge states are 
generated. 

For < /i < Ai/2, edge currents arc propagating 
in opposite directions for opposite spin polarizations at 
both sample boundaries [sec Fig.[5Jal)], which is similar 
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FIG. 3: (Color online) Similar to that of Figs. O but for four representative configurations of two interfaces (a) Case III with 
Bi = 2 T, B2 = 3 T and Bs = 4 T, (b) Case IV with Bi = 3 T, B2 = 2 T and B3 = 4 T, (c) Case V with Bi = 2 T, 
Ba = -4 T and B3 = 3 T, and (d) Case VI with Bi = 2 T, B2 = 4 T and B3 = -3 T. 



to the quantum spin Hall state in graphene by consid- 
ering intrinsic SO interaction^ and also similar to the 
state predicted in Refl23 in a uniform magnetic field. For 
fi > A2/2, however, the edge currents at sample bound- 
aries become normal because they propagate in the same 
direction for both spin polarizations [see Fig. [5][a3)]. In- 
terestingly, when the chemical potential fi lies within the 
gap Ai/2 < /i < A2/2, the interface channels are opened 
and two spin-down edge states with the same flowing di- 
rection emerge at the interface, and certainly, they are 
FSP [Fig. [2Ja2)]. Meanwhile, the edge current at one 
sample boundary becomes normal whereas that at the 
other boundary is still left novel. For case where the 
magnetic fields in the two sides of the interface are in 
opposite directions, similar conclusion can be made ex- 
cept that when fj, > A2/2, no FSP state is present but 
the interface current channel is not closed and there ex- 
ist a normal edge current at the interface instead [see 
Figs.[lbl)-lli;b3)]. 

From the viewpoint of application, one can use the 
properties mentioned above to manipulate the FSP edge 
currents. Imaging placing a voltage gate on the top of 



graphene sample, and by changing gradually the voltage, 
one can vary the carrier density in the sample, so the 
chemical potential fi gradually. We take Case I as an 
example. When fi lies within the minimum spin gap \fJ.\ < 
Ai/2, the interface edge current channel is closed. As 
increases, the interface channel with a FSP edge current 
is opened if fi lies within the two gaps Ai/2 < < 
A2/2. Finally the channel is closed again in the interval 
/i| > A2/2. This implies that one can control the FSP 
edge currents electrically. 

The spin gap due to Zeeman splitting is estimated as 
3 K at B = 2 T, whereas the corresponding first 
LL El is estimated to be about 500K. Furthermore, it 
was argued that the spin gap can be enhanced to a few 
hundred kclvin for a realistic magnetic field due to elec- 
tron exchange interaction)^ So the FSP current is ex- 
pected to be observed experimentally in relatively large 
range of temperatures. To visualize the effect of the 
Zeeman splitting clearly and not changing the physics, 
relatively large spin gaps have been chosen for calculat- 
ing the energy spectrum in Fig. [5] We also note that 
the localization character along x direction for the edge 



4 



states i s well d efined beeause the magnetic length here is 
Ib = ^hc/eB = 18 nni, much less than the width of the 
sample L (=213 nm). 

We next consider the case of NMF with two interfaces. 
There are twelve in total different configurations for pat- 
tern of magnetic field, but among them, only four( Case 
III- Case IV) are representative and of great interest, 
which is shown in Fig. [3] Three spin gaps are formed 
in each case, so similar to Fig. [2l the states between the 
zero LLs can be generally classified into four regimes, 
according to the position of the chemical potential /x. 

Case III is very special. For this interesting case, if 
generally, there are more interfaces with the magnetic 
field in each region in the same direction and increasing 
monotonically, the FSP interface edge channels will be 
opened one by one from the low-field interface to high- 
field one, where the old one closing and new one opening 
happen simultaneously. For Case IV, when /i increases 
gradually from < /i < Ai/2 to /z > A3/2, the two 
interface channels are opened simultaneously, after that 
one channel is closed and finally the other one is closed. 
For Case V and Case VI, although the FSP interface 
channels are opened and closed in different orders, these 
channels are not closed in the end but become normal 
ones instead [Figs. [31; c4) and 121 d4)]. 

Based on the above discussion, we come to a general 
conclusion that for a NMF with many interfaces, so long 
as the magnetic length is much less than the width of 
each uniform region, there always exist FSP edge states 
either at all the interfaces or at several of them, when the 
chemical potential ji lies within the gap Ami„/2 < \^\ < 
^max/'2- AH the calculations above are performed for 
the zigzag graphene ribbon. For the armchair graphene 
ribbon, similar edge states and so similar conclusion can 
be obtained. 



B. Rashba spin-orbit coupling effect 
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FIG. 4: (Color online) (a) Electron energy spectrum of 
graphene in the presence of Rashba SO interaction for Case I. 
Parameters are the same as that in Fig. [2] (a) and A = Ai/50. 
(b) The spin tilt angle as a function of E for different edge 
states. The inset in (b) is an enlarged view of the tilt an- 
gle within the spin gap — A2/2 < fi < — Ai/2. Open and 
filled triangles denote the interface edge states, while open 
and filled circles denote the sample boundary edge states. 



continuously with ky. Therefore, theoretically speaking, 
the FSP interface current is still preserved at T = K 
[Fig. m (a)], except that the edge current is now polar- 
ized along other direction instead which slight deviates 
from z axis by an angle 9{kp). In the spin gap region 
Ai/2 < \fi\ < A2/2, for the parameters shown in Fig. 4, 
the tilting angle for the interface edge states is estimated 
to be |6'| < 0.0074. ft is very small and changes nearly 
linearly with energy with the slope A6/AE ~ 1.3 in this 
gap [see the inset in Fig. |3] (b)]. Physically, the slope 
A9/ AE is believed to be nearly proportional to A, which 
is rather small in real situations (For a perpendicular 
electric field E, - 50V/300 nm^i22 A ~ 0.07 K,^^). At 
finite temperatures, in order to obtain a almost FSP in- 
terface current, the thermodynamical energy ksT should 
be much less than the spin gap (A2 — Ai)/2, which is es- 
timated to be about 30 K for B = 2 TM Therefore, in 
many situations achievable in experiments, the spin-flip 
processes induced by the Rashba SO interaction can be 
negligible, and an almost FSP current can be realized. 



We now discuss the effect of the spin-fiip processes 
caused by spin-orbit(SO) interaction. To estimate the 
influence of the SO interaction on the polarized edge cur- 
rent, we only consider here the Rashba SO interaction as 
an example, whose Hamiltonian can be written as 

Hr^xY, W ^ ^■'"c| (a X d,, ) • zcj + H.c] , (4) 

<i]> 

where A is the Rashba SO coupling strength and a is the 
Pauli matrix, d^^ is a vector pointing from j to its nearest 
neighbor i, and z is a unit vector in z direction. Figure 
|4] (a) shows the spectrum at low energies with Rashba 
SO interaction, which is found to be almost unchanged 
compared with the case without SO coupling [Fig. [2] (a)]. 
However, in contrast to the case without SO interaction 
where electrons are polarized along z axis, the spin po- 
larization of the edge states in the presence of the SO 
interaction can be described by S = (sin0, 0, cos0) with 
the tilting angle deviated from the z axis, changing 



C. Exporting the FSP edge currents 

So far, we have shown explicitly how to generate and 
manipulate the FSP edge currents, but another associ- 
ated and interesting question which naturally arises is 
still unsolved, which is how to export the FSP current 
to serve as a current source for a spin related device. 
This is not a trivial question, because to export the cur- 
rent, one has to destroy the interface in some extent 
to introduce another lead, but the FSP current is com- 
posed of edge states, which will disappear when the in- 
terface is destroyed or removed away. For this purpose, 
we design an experimental setup shown schematically In 
Fig. [S{a), which is expected to export FSP current to 
the possible spin related device connected with the arm- 
chair graphene sample by a narrow graphene strip. Here 
the armchair graphene ribbon is adopted, because for 
the zigzag graphene ribbon, the existence of the flat edge 
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(a) (c) 




FIG. 5: (Color online) (a) Schematic diagram of our proposed 
experimental setup for application in spintronics device, (b) 
The corresponding eigenvalues calculated on the sample in (a) 
with the magnetic flux per hexagon (j) = 2-k and 1 . 5 x 27r/80 
(B ~ 1000 T). The sample size is 40 x 120 consisting of 40 
armchair chains with 120 sites in each chain, while the narrow 
graphene strip size is 4 x 40. The Zeeman energy is taken as 
hz = 0.1 and hz ~ 0.15 for the two regions respectively, (c) 
The corresponding three-terminal device with the sample size 
as (b). (d) Geometry of the graphene sample connected with 
three leads for calculating transmission coefficients. To reduce 
the scattering between the sample and the leads, there exist 
a transition region with length 1.14 nm, where the magnetic 
field is made to be linearly decreased to zero when approach- 
ing leads. 

bands will cause ferromagnetic instability and lead to lo- 
cal ferromagnetic polarizations at sample edgeS)^ giving 
rise to spin-flip processes. 

To testify this idea, we first calculate the eigenvalues 
of this particular system for an open boundary condi- 
tion. We find that LLs arc still present even in the 
condition that the interface is partially destroyed by the 
introduced lead. However, states for both spin species 
exist between the LLs Ai/2 < |^| < A2/2 [Fig. [DJb)], 
leaving the question whether spin-up electrons in this 
regime contribute to the current flowing through the de- 
vice still unaddressed. Therefore, in order to conflrm 
that the current flowing through the device is FSP, we re- 
fer to the three-terminal setup as illustrated in Fig. [Sfc) , 
where the rectangular graphene sample is connected with 
three ideal semi-infinite leads. The detail of these con- 
nections between graphene sample and the leads is shown 
in Fig. [5jd). In this situation, four bonds across the 
interface are broken, so that lead 3 made of the nar- 
row graphene strip can be connected with the sample. 
This strip is then upright so as to avoid the influence 
of magnetic field. According to the Fisher-Lee relation, 
the electron transmission coefficient from lead q to lead 
p for multi-terminal Landauer-Biittiker formula is given 

=Tr[r^G«r^G^], (5) 

The scattering rate matrix Fp due to the coupling to lead 
p has a form 

r^(z,j) = z[I]^,(*,j)-S^<.(*,j)], (6) 



where = [S^^]''' is the self-energy due to lead p and 
the elements of the matrix 'S^^ are given by 



with the function x„i(«) = yj^ifSm^^ describing 
the transverse profile of mode m in lead p and the 
wave vector determined by the equation p, = —2 cos k,n — 
2 cos ■ Here TV is the number of sites connected with 
the leads and m = 1,2,-- - ,A^. The retarded Green's 
function of the graphene sample is written as 

G« = [i?/-i7o-S«]-\ S^=^E^, (8) 

p,(j 

and ~ [G^]^. Here, to make the calculation simple, 
the three leads are assumed to be described by a nearest- 
neighbor tight-binding model on a square lattice with the 
same hopping integral as graphene. 

The transmission coefficients T^^ are plotted in Fig. [5l 

For spin-up electrons, T^2 ^.nd show quasi-periodic 
oscillation and the maxima of them are the quantized 
value 1 . This oscillating behavior is similar to that found 
in graphene-based ferromagnetic double junctions^and 
can be attributed to electrons' multiple reflections be- 
tween terminal 1 and 2. As the energy increases gradually 
from to 0.2, both ^i^d are approaching to 1, sug- 
gesting one left-moving spin-up edge state on one sample 
boundary and one right-moving one on the other, sharing 
the similar picture with Fig.[2Ia). However, most impor- 
tantly, we find that within the error range, T^g = T^-^ ~ 
and tJj ~ — in the whole energy range < i? < 0.2 
we are interested in. This is corresponding to the absence 
of spin-up edge states at the interface. 

For spin-down electrons, when < yU < Ai/2 and 
/i > A2/2, ~ T21, the situation is similar to that 
for spin-up electrons. However, when Ai/2 < p < A2/2, 
T^2 ^-nd T21 show strongly fluctuations, and deviate from 
the quantized value 2 [Fig.tHl^a)], which should be taken 
in the ideal case shown in Fig.[2Ja). This may be ascribed 
to the bonds breaking at the interface of terminal 3, and 
lattice mismatching at the interfaces of terminal land 2 
which results in electrons scattering from edges states to 
the nearby zeroth LL [see Fig.OJb)]. In the same range, 

the fact that T^^ > Tg^ > and T^2 > ^^3 > indi- 
cates that both the left- and right-moving currents exist, 
but the net current is left moving. The finitencss of 
also implies that a small fraction of spin-down electrons 
incoming from terminal 3 have to be scattered to sam- 
ple boundaries via bulk states, since there is no current 
channel at the interface from terminal 3 to 2. Similarly, 
the finiteness of T^^ > means that the right-moving 
spin-down electrons propagating along sample bound- 
aries have some probability to be scattered to terminal 3 
via bulk states. These processes arc schematically shown 
in Fig. El^g). We remark that although these scattering 
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FIG. 6: (Color online) (a)-(f)The transmission coefficients as 
functions of the Fermi energy. The coefficients in (a), (c) and 
(e) are for spin-down electrons, while the others for spin-up 
electrons. The dash lines at S = 0.1 and 0.15 indicate the 
positions of two spin gaps, (g) Schematic diagram of electron 
current flow when Ai/2 < < A2/2 (i.e., 0.1 < /i < 0.15). 
(h) The current following out of terminal 3, I3, when terminal 
voltages are chosen as Vi = V2 = V and V3 = 0. The inset 
in (a) is an enlarged view of T^2 ^21 S^P < < 

Ai/2, and the inset in (b) is an enlarged view of Tjt, and tJ^ 
in the energy range < iJ < 0.2. 



processes will cause resistance, spin is conserved since 
there is no spin-flip processes considered here. This leads 
to the conclusion that the electron current flowing out of 
terminal 3 has spin-down polarization, and so is FSP. 

In actual application, one can apply the electrical volt- 
age to the three terminals as Vi = V2 = ^ and V3 — 0, 
then the FSP (spin-down) electrons will be transported 
from lead 3 to leads 1 and 2. In this situation, the FSP 
current flowing out of terminal 3, I3 ~ {T^^ + T^2)^^'^ 
is estimated as OMVe'^/h - 0.95Ve'^/h [see Fig. i^h)]. 
Indeed, the width of the strip is d = 0.87 nm and the 
magnetic length is about Ib ~ 0.8 nm. The fact of d ^ Ib 
guarantees that there are sufficient electrons through the 
strip. Therefore, the graphene sample with the designed 
setup under a NMF can be viewed as the FSP current 
source for a spin related nanodevice. 

IV. SUMMARY 

In summary, we have numerically investigated the FSP 
edge currents at the interface in graphene under a nonuni- 
form magnetic field. The manipulation of the FSP cur- 
rents is realized by a gate voltage, and its potential ap- 
plication in spintronic devices is also discussed. An ex- 
perimental setup is designed and is expected to export 
the FSP current, which is confirmed numerically by a 
three-terminal transmission calculation. 



Acknowledgments 

The author J. An thanks X. G. Wan, H. H. Wen, J. 
X. Li, Q. H. Wang for useful discussions. This work 
was supported by NSFC Projects No. 10504009, and 
No. 10874073, and 973 Projects No. 2006CB921802, No. 
2006CB601002, and No. 2009CB929504. 



* Electronic address: anjin@nju.edu. cn| 

t Electronic address: cdgongsc@nju.edu.cn' 

^ M. Zarea and N. Sandler, Phys. Rev. B 79, 165442 (2009). 

^ Y. K. Kato, R. C. Myers, A. C. Gossard, D. D. Awschalom, 

Science 306, 1910 (2004). 
^ J. Sinova, D. Culcer, Q. Niu, N. A. Sinitsyn, T. Jungwirth, 

and A. H. MacDonald, Phys. Rev. Lett. 92, 126603 (2004). 

* C. L. Kane and E. J. Mele, Phys. Rev. Lett. 95, 226801 
(2005). 

^ R. A. de Groot, F. M. Mueller, P. G. van Engen, and K. 

H. J. Buschow, Phys. Rev. Lett. 50, 2024 (1983). 
® J-H. Park, Nature 392, 794 (1998). 

K. R. Kittilstved, D. A. Schwartz, A. C. Tuan, S. M. Heald, 

S. A. Chambers, and D. R. Gamelin, Phys. Rev. Lett. 97, 

037203 (2006). 

® A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. 
Novoselov, and A. K. Geim, Rev. Mod. Phys. 81, 109 
(2009). 



^ A. K. Geim and K. S. Novoselov, Nat. Mater. 6, 183 (2007). 
1° J. C. Meyer et al.. Nature (London) 446, 60 (2007). 
" S. Cho, Yung-Fu Chen, and M. S. Fuhrerb, Appl. Phys. 
Lett. 91, 123105 (2007). 

N. Tombros, C. Jozsa, M. Popinciuc, H. T. Jonkman, and 

B. J. van Wees, Nature (London) 448, 571 (2007). 
1^ E. W. HiU, A. K. Geim, K. S. Novoselov, F. Schedin, and 

P. Blake, IEEE Trans. Magn. 42, 2694 (2007). 
" M. B. Lundeberg and J. A. Folk, Nat. Phys. 5, 894 (2009). 

Y-W. Son, M. L. Cohen, and S. G. Louie, Nature. 444, 

347 (2006). 

J. Jung, T. Pereg-Barnea, and A. H. MacDonald, Phys. 
Rev. Lett. 102, 227205 (2009). 
" J. Jung and A. H. MacDonald, Phys. Rev. B 79, 235433 

(2009) . 

L. Xu, J. An, and C.-D. Gong, Phys. Rev. B 82, 155421 

(2010) . 

M. Ramezani Masir, P. Vasilopoulos, A. Matulis, and F. 



7 



M. Peeters, Phys. Rev. B 77, 235443 (2008). 
^° S. Park and H.-S. Sim, Phys. Rev. B 77, 075433 (2008). 
M. Cerchez, S. Hugger, T. Heinzel, and N. Schulz, Phys. 
Rev. B 75, 035341 (2007), and references therein. 
M. L. Leadbeater, C. L. Foden, J. H. Burroughes, M. Pep- 
per, T. M. Burke, L. L. Wang, M. P. Grimshaw, and D. A. 
Ritchie, Phys. Rev. B 52, R8629 (1995). 
S. J. Bending, K. von Klitzing, and K. Ploog, Phys. Rev. 
Lett. 65, 1060 (1990). 

D. A. Abanin, P. A. Lee, and L. S. Levitov, Phys. Rev. 
Lett. 96, 176803 (2006). 
■^^ D. Huertas-Hernando, F. Guinea, and A. Brataas, Phys. 
Rev. B 74, 155426 (2006). 

H. Min, J. E. Hill, N. A. Sinitsyn, B. R. Sahu, Leonard 
Kleinman, and A. H. MacDonald, Phys. Rev. B 74, 165310 



(2006). 

K. S. Novoselov, A. K. Geim, S. V. Morozov, D. Jiang, Y. 
Zhang, S. V. Dubonos, I. V. Gregorieva, and A. A. Firsov, 
Science 306, 666 (2004). 

M. Fujita, K. Wakabayashi, K. Nakada, and K. Kusakabe, 
J. Phys. Soc. Jpn. 65, 1920 (1996). 

L. Sheng, D. N. Sheng, and C. S. Ting, Phys. Rev. Lett. 
94, 016602 (2005). 

E. M. Hankiewicz, L.W. Molenkamp, T. Jungwirth, and J. 
Sinova, Phys. Rev. B 70, 241301(R) (2004). 
C. Bai, J. Wang, S. Jia, and Y. Yang, Appl. Phys. Lett. 
96, 223102 (2010). 

E. Prada, P. San- Jose, and L. Brey, Phys. Rev. Lett. 105, 
106802 (2010). 



